Logarithmic rate dependence in deforming granular materials 
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Rate-independence for stresses within a granular material is a basic tenet of many models for 
slow dense granular flows|l|, ^l_^i.M_3- contrast, logarithmic rate dependence of stresses is 
found in solid-on-solid frictioii|^^g|7 in geological settings[|, 0, and elsewhere [|ll|, ^ |l^, 
In this work, we show that logarithmic rate-dependence occurs in granular materials for plastic 
(irreversible) deformations that occur during shearing but not for elastic (reversible) deformations, 
such as those that occur under moderate repetitive compression. Increasing the shearing rate, Q, 
leads to an increase in the stress and the stress fluctuations that at least qualitatively resemble 
what occurs due to an increase in the density. Increases in Q also lead to qualitative changes in 
the distributions of stress build-up and relaxation events. If shearing is stopped at t = , stress 
relaxations occur with o{t)/a[t = 0) ~ A\og{t/to) ■ This collective relaxation of the stress network 
over logarithmically long times provides a mechanism for rate-dependent strengthening. 

PACS numbers: PACS numbers: 46.10.+Z, 47.20.-k 



Slow granular flows, the subject of this letter, are typ- 
ically described in the context of Mohr- Coulomb friction 
models Q that resemble those used for describing friction 
between two solid bodies [T5[| . In the well-known solid 
friction scenario an object on a frictional sur- 
face will resist a force and remain at rest provided the 
magnitude of the tangential force is less than the prod- 
uct of a static friction coefficient and the normal force. 
This picture was translated into the granular context (for 
dense gra nular systems characterized by networks of force 
chains [p^ ) by Coulomb ||l9|| and more recent authors [0, ||, 
where the normal and tangential forces are replaced by 
corresponding normal and shear stresses, and the surface 
of interaction is replaced by a plane within the material. 
For large enough tangential force (shear stress) relative to 
the normal force (normal stress) sliding friction (failure 
in a granular material) occurs. In these pictures, sliding 
friction (deformation following failure) is independent of 
the speed of sliding (the shear rate). 

In reality, experiments in diverse contexts have shown 
that solid friction exhibits a logarithmic dependence on 
rate and that static frictional contacts strengthen log- 
arithmically with age. These experiments span a vast 
range of lengths, and include studies at the atomic JT^, 
lab^, I and geological scale |, |lO). Recent 

experiments by Ovarlez et al. ^] using granular materials 
sliding against the interior wall of a piston showed clear 
rate dependence that was associated with aging effects of 
individual solid-friction contacts and with the force net- 
work. Additionally, experiments by Nasuno et al. j^, ^ in 
which a solid surface was pushed across a granular bed 
showed a slow strengthening (or aging) with time of the 
(quasi-static) force. In the present experiments, which 
zoom in uniquely on the grains, we show that there is a 
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logarithmic rate-dependence in slowly sheared granular 
materials that is associated with irreversible rearrange- 
ments of the grain contacts within the material itself. 
This is manifest as a strengthening with rate, and it is 
unique to granular systems and possibly other jammed 
sy stems [pO|; i.e., it does not depend per se on the loga- 
rithmic strengthening of frictional contacts. 

The experiments described here were carried out with 
a 2D realisation of a granular system. The particles were 
made of a photoelastic material and were either relatively 
thin disks or flat particles with a pentagonal cross section. 
By using a photoelastic material, we could determine 
forces at the grain scale. More detailed descriptions of the 
experimental apparatus and methods have been given 
elsewhere for related experiments , and here we 

provide only essential information for two qualitatively 
different experiments. For each experiment, the humid- 
ity varied by no more than ±5%. In the first of these 
two experiments, the photoelastic particles were sheared 
in an annular geometry (inset. Fig. |^) that was bounded 
on the inside by a rough wheel, the source of shearing, 
and on the outside by an equally rough static ring. The 
particles sat on a horizontal, flat powder-lubricated Per- 
spex sheet, and were contained from above by a similar 
sheet. The whole was placed in a polariscope consisting 
of a pair of right and left circular polarizers, a light source 
and a video camera. Photoelastic images were obtained 
at rates up to 15 Hz, captured by a framegrabber, and 
processed on-the-fly to determine the forces on the par- 
ticles within the field of view of the camera (~ 100 to 
^ 200 particles). Specifically, we computed the force on 
each particle in each frame and summed over all parti- 
cles in a frame. We stored this integrated force for each 
image to produce a time series over very long sampling 
times, between 2^^ and 2^^ points depending on the rate 
of shearing. To a reasonable approximation (^ 10%), this 
quantity is proportional to the pressure, and for simplic- 
ity, we will refer to it as the stress. 
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FIG. 1; Time series for stress (in arbitrary units) for a range 
of shearing rates Q that span the experimental range. These 
data are for pentagonal particles, but data for disks are qual- 
itatively similar. 



Key control parameters are the rate of shearing, Q — 
2tt/T, where T is the period for one rotation of the shear- 
ing wheel, and the density of the system, which we give 
in terms of the area fraction, 7, occupied by the par- 
ticles. As shown previously 1^ |2|], there is a critical 
value of 7c, below which the shearing ceases, and we typ- 
ically reference 7 to this value. (7c depends on prop- 
erties such as the particle shape and possibly the ge- 
ometry of the container.) It is also useful to give the 
rotation rate, / = 1/T — U/2tt. In these experiments, 
0.03 mHz < / < 60 mHz. 

In Fig. |l|, we show several examples of stress time series 
for n's ranging over about three orders of magnitude in 
the shearing rate. The three images of this figure give an 
indication of the effect of shearing rate: the time series 
for the slowest rate is clearly of lower amplitude and more 
intermittent than that for the highest rate. 

The rate dependence is seen clearly in the mean (time- 
averaged) stress vs. fi. Fig. ^. These data are consistent 
with a logarithmic rate dependence for all 7's. Similar ef- 
fects have been seen in solid-on-solid friction experiments 
but with a significant difference: in the present experi- 
ments, there is a logarithmic strengthening of the mean 
stress with rate, while other experiments [|[ ||, |ll| show a 
logarithmic strengthening with waiting time and/or hu- 
midity. If we fit these data to the form a = A\og{il./ilo), 
then the amplitude. A, is sensitive to the density, but the 
characteristic frequency, fio, depends only weakly on 7, 
and corresponds to a typical time of 27r/f2o w 30 days. 
However, there is no compelling reason to believe that 
this relation holds in the limit ^ 0. A more realis- 
tic description may he a = Alog[{fl + f7i)/$lo], where 
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FIG. 2: Mean stress vs. shearing rate, Q for a range of den- 
sities for shearing. The lin-log scales emphasize the fact that 
the data is consistent with a logarithmic variation of the stress 
with Q. Inset: sketch of Couette shear apparatus. 



If a sample is sheared under steady state conditions 
and the shearing wheel is abruptly stopped, the stress 
relaxes over very long time scales. In Fig. y[ we show 
the relaxation of stress vs. time on a semi-log plot. 
Note that there are two time scales: a quick initial decay 
lasting 20 s and a much longer process such that the 
stress network may still be relaxing after 20 hrs. Fig. ^ 
shows fits to the stress decrease with time, of the form 
a{t = 0)/a{t) ~ A\og{t/to). Although the relaxation is 
not completely uniform in time, a logarithmic functional 
form is not unreasonable for many cases (after an initial 
rapid relaxation). This was not always the case, and an 
interesting example of a major relaxation event is seen 
in one of these runs about 350 s after the stopping of 
shearing; thereafter, continued slow relaxation occurred. 
A key insight is that, while the micro-contacts between 
particles may be strengthening in time|||, the 
stress network as a whole is relaxing over very long time 
scales. 

We gain additional insight into the effect of rate 
changes by determining the distribution of build-up 
events (monotonic increases in stress) and release events 
("avalanches" - monotonic decreases in stress). Here, 
there are two aspects that are of interest: the size of the 
build ups and avalanches. Act, and the angle, 9, through 
which the shearing wheel turns during the course of a 
build-up or an avalanche. In the top part of Fig. ^, we 
show an example of an "avalanche" event, and in the bot- 
tom, we show distributions for Act and 9 for avalanches. 
Note that for near-critical 7's and for the slowest fi's, the 
avalanche distributions are consistent with power laws. 
At higher 7 and f2 the distributions for both build-up (not 
shown due to space constraints) and avalanche events are 
roughly exponentials. Interestingly, both increases in 
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FIG. 3: Relaxation of stress vs. time of static configura- 
tions of grains, previously undergoing plastic deformation due 
to shear. The dashed lines are a fit to a{t)/a(t = 0) ~ 
A\og{t/to) 
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FIG. 4: Top: a short segment of a time series that identifies 
a single avalanche, Act and 9. Bottom: distributions of stress 
avalanches for a dense system of pentagons, 7 = 7c + 0.0226, 
for several rates on a log-lin graph (left) and a log-log graph 
(right). 



and in 7 tend to have the same quahtative effect on the 
statistics of avalanches and build-ups. 

In order to further explore the origin of this rate depen- 
dence, we have carried out a second set of experiments, 
using the same photoelastic techniques, in which the par- 
ticles were confined to a box with three fixed sides. The 
sample was then compressed and released repeatedly on 
the fourth side by means of an oscillating piston, as in 
the sketch of the inset of Fig. |[ The driving velocity of 
the piston was chosen to match the corresponding veloc- 



FIG. 5: Mean stress vs. oscillation rate, SI for several densi- 
ties for oscillatory compression. Data for both disks and pen- 
tagons is shown. The lin-log scales emphasize the fact that 
there is no variation of the stress with Q within experimental 
uncertainties. Inset: sketch of piston apparatus. 



ities of the shear wheel in the first experiment. The key 
difference between the two sets of experiments was that 
in the first set there was continuous plastic deformation 
while in the latter set there was none; i.e., the particles 
maintained their relative position. As shown in the body 
of the figure, which gives a vs. fi, there was no rate 
dependence of this process, within experimental error. 

To conclude, we have demonstrated that sheared gran- 
ular materials exhibit logarithmic rate dependence that is 
tied in an essential way to particle rearrangements (plas- 
tic deformation). The nature of this rate dependence is 
further illuminated by the slow relaxation that occurs 
when shearing is halted. The latter experiments in par- 
ticular suggest that there are slow collective rearrange- 
ments of the particles that can occur over large time 
scales. This slow relaxation of stresses over long times 
is consistent with increases in with rate: force chains 
generated by shearing cannot completely relax on the 
time scales over which new chains are formed, an effect 
that is exacerbated by increasing fl. It seems likely that 
this effect is intimately related to the long time scales 
reported for compaction in granular materials |Q. For 
compaction, collective rearrangements become progres- 
sively more difficult over time because each new rear- 
rangement requires the involvement of larger collections 
of particles. Similarly, in the present experiments, re- 
laxation events involving collections of particles become 
less probable over time due to geometric effects and to 
the fact that the available elastic energy is gradually re- 
duced with each rearrangement. An interesting question 
is whether similar properties are seen in other jammed 
systems [pOl such as colloids or foams. 



4 



Acknowledgments als Science & by NASA. We appreciate helpful conversa- 

tions with Prof. Sue Coppersmith. 

Acknowlegements: This work was supported by NSF 
Division of Materials Research, NSF Division of Materi- 



[1] Wood, D. M. Soil Behaviour and Critical State Soil 
Mechamcs (Cambridge University, Cambridge, England, 
1990). 

[2] Nedderman, R. M. Statics and Kinematics of Granular 
Materials (Cambridge Univ. Press, Cambridge, 1992). 

[3] Jaeger, H. M., Nagel, S. R. & Behringer, R. P. Granular 
solids, liquids, and gases. Rev. Mod. Phys. 68, 1259-1273 
(1996). 

[4] Herrmann, H. J., Hovi, J. P. & Luding, S., eds. Physics of 
Dry Granular Media, (NATO ASI Series, Kluwer, 1997). 

[5] Behringer, R. P. & Jenkins, J. T., eds. Powders and 
Grains 97, (Balkema, Rotterdam, 1997) 

[6] Ovarlez, C, Kolb, E. & Clment, E. Rheology of a 
confined granular material. Phys. Rev. E 64, 060302- 
060305(R) (2001). 

[7] Nasuno, S., KudroUi, A., Bak, A. & Gollub, J. P. Time- 
resolved studies of stick-slip friction in sheared granular 
layers. Phys. Rev. E 58, 2161-2171 (1998). 

[8] Losert, W., Geminard, J.-C, Nasuno, S. & Gollub, J. P. 
Mechanisms for slow strengthening in granular materials. 
Phys. Rev. E Ql, 4060-4068 (2000). 

[9] Dieterich, J. H. Modeling of rock friction 1: Experimental 
results and constitutive equations. J. Geophys. Res. 84, 
2161-2168 (1979). 
[10] Ruina, A. L. Slip instability and state variable friction 

laws. J. Geophys. Res. 88, 10359-10370 (1983). 
[11] Berthoud, P., Baumberger, T., G'Sell, C. & Hiver, J. M. 
Physical analysis of the state- and rate-dependent friction 
law: Static friction. Phys. Rev. B 59, 14313-14327 (1999). 
[12] Baumberger, T., Berthoud, P. & Caroli, C. Physical anal- 
ysis of the state- and rate-dependent friction law. II. Dy- 
namic friction. Phys. Rev. B 60, 3928-3939 (1999). 



[13] Carlson, J. M. & Batista, A. A. Constitutive relation for 
the friction between lubricated surfaces. Phys. Rev. E53, 
4153-4165 (1996). 

[14] Gnecco, E. et al. Velocity dependence of atomic friction. 
Phys. Rev. Lett. 84, 1172-1175 (2000). 

[15] Heslot, F., Baumberger, T., Perrin, B., Caroli, B. & Car- 
oli, C. Creep, stick-slip, and dry-friction dynamics: Ex- 
periments and a heuristic model. Phys. Rev. E 49, 4973- 
4988 (1994). 

[16] Persson, B. N. J. Sliding Friction: Physical Principles 
and Applications (Springer- Verlag, Berlin, 1998). 

[17] Meyer, E., Overney, R. M., Dransfeld, K. & Gyalog, T. 
Nanoscience, Friction and Rheology on the Nanometer 
Scale (World Scientific, Singapore, 1998). 

[18] Howell, D., Vcje, C. & Behringer, R. P. Stress fluctua- 
tions in a 2D granular Couette experiment: A continuous 
transition. Phys. Rev. Lett. 82, 5241-5244 (1999). 

[19] Coulomb, A. Memoires de Mathematiques et de Physique, 
(L'Imprimerie Royalc, Paris, 1776). 

[20] Liu, A. L. & Nagel, S. R. Nonlinear Dynamics: Jamming 
is just not cool anymore. Nature. 396, 21-22 (1998). 

[21] Heywood, R. B. Designing by Photoelasticity (Chapman 
and Hall Ltd., London, 1952). 

[22] Howell, D. W., Behringer, R. P. & Veje, C. T. Fluctua- 
tions in granular media. Chaos 9, 559-572 (1999). 

[23] Gcng, J. et al. Footprints in sand: The response of a 
granular material to local perturbations. Phys. Rev. Lett. 
87, 035506-035509 (2001). 

[24] Nowak, E. R., Knight, J. B., Ben-Naim, E., Jaeger, H. M. 
& Nagel, S. R. Density fluctuations in vibrated granular. 
Phys. Rev. E57, 1971-1982 (1998). 



